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O ; Abstract 

o 

: The production rates of D*±, D*±, D±, D7D0, D^^, and A+M^" in Z ^ cc 

I decays are measured using the LEP I data sample recorded by the ALEPH detector. 

Q ' The fractional energy spectrum of the D*^ is well described as the sum of three 

p • contributions: charm hadronisation, b hadron decays and gluon splitting into a pair 

1^ ■ of heavy quarks. The probability for a c quark to hadronise into a D*"*" is found to be 
/(c ^ D*+) = 0.233 ± 0.010(stat.) ± 0.011(syst.). The average fraction of the beam 
energy carried by D*^ mesons in Z cc events is measured to be {Xe{^*^))^^ = 

X ; 0.4878 ± 0.0046(stat.) ± 0.0061(syst.). The D*± energy and the hemisphere mass 

' imbalance distributions are simultaneously used to measure the fraction of hadronic Z 



decays in which a gluon splits to a cc pair: Ug^cc = (3.23 ± 0.48(stat.) ± 0.53(syst.))%. 
The ratio of the Vector/ (Vector+Pseudoscalar) production rates in charmed mesons is 
found to be Py = 0.595 ± 0.045. The fractional decay width of the Z into cc pairs is 
determined from the sum of the production rates for various weakly decaying charmed 
states to be Rc = 0.1738 ± 0.0047(stat.) ± 0.0116(syst.). 
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1 Introduction 



Charm quarks are produced in 40% of the hadronic Z decays, mostly coming from 
Z ^ cc or bb, but also from gluon splitting to cc or bb. This paper presents a wide- 
ranging survey of the production of the weakly decaying charm states D°, D+, D+ 
mesons and baryons and the first excited states D*+, D*+ in Z decays (here and 
throughout this paper, charge conjugation is implied). It extends the work of previous 
publications to the full data sample collected by ALEPH between 1991 and 1995, 
which consists of four million hadronic Z decays. A new algorithm based on the mass 
of particles from secondary vertices greatly improves the separation of charm states 
resulting from bb and cc production. 

The paper is organised as follows. After a brief description of relevant details of 
the ALEPH detector in Section 2, the procedures adopted to reconstruct the above 
charm states are given in Section 3. The cleanest charm signal is obtained for the D*"*" 
mesons and Section 4 is devoted to the fractional energy spectrum for their production, 
which now covers the whole range from 0.1 to 1.0, and the relative strengths of the 
bb, cc and gluon splitting components. In Section 5 a new technique involving the 
heavy and light hemisphere masses is introduced in order to obtain a more accurate 
measure of the rate for gluon splitting to charm quarks. A measurement of the D*^ 
production rate is given in Section 6. Determinations of the production rates for the 
weakly decaying charm states are presented in Section 7 and this information is used 
in Section 8 to yield a measurement of Rc = r(Z cc)/r(Z hadrons) from charm 
counting. 

2 The ALEPH detector 

A detailed description of the ALEPH detector and its performance can be found in 
Refs. 0. Only a brief review is given here. 

Charged particles are detected in the central part, consisting of a two-layer silicon 
vertex detector with double-sided (r-0 and z) readout, a cylindrical drift chamber and 
a large time projection chamber (tpc), which together measure up to 33 coordinates 
along the charged particle trajectories. Tracking is performed in a 1.5 T magnetic field 
provided by a superconducting solenoid. For high momentum tracks the combined 
system yields a 1/pt resolution of 6 x 10~^(GeV/c)^^ and an impact parameter 
resolution of 25 /xm in both the r-0 and z projections. The TPC also provides up to 338 
measurements of ionization [dE/dx) allowing particle identification to be performed. 
For tracks having at least 50 energy deposition measurements, the variable Xh is defined 
as {Imeas — Ih)/o'h, where /meas is the measured value of dE/dx, Ih is the expected value 
for particle type h [h = p,K or tt), and ah is the expected uncertainty. 

Data from the vertex detector is particularly important for b tagging for which a 
lifetime-mass algorithm has been designed to provide discrimination between b quarks 
and lighter quarks. This algorithm uses the significance S of the three-dimensional 
impact parameter for each charged track to define the confidence level for it to be 
consistent with coming from the primary vertex. These confidence levels are combined 
to obtain the confidence level V-h for all tracks of an hemisphere to come from the 
primary vertex. The distribution of V-h is strongly peaked near zero for b hemispheres. 
A second b tagging variable makes use of the b/c hadron mass difference: the tracks 
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in a hemisphere are grouped together in order of decreasing S until an effective mass 
of 1.8 GeV/c^ is reached for the system, and the confidence level fi-j^ for the last 
added track to come from the primary vertex is used as the tagging variable. This 
is peaked at zero for b quarks. These two tagging variables are combined in a single 
one: Btag = —{0.7logiofiT-i + O.SlogioV-H- A detailed description of this algorithm can 
be found in Ref. [^. 

The electromagnetic calorimeter is a lead/wire- chamber sandwich operated in 
proportional mode. It is read out in projective towers of typically 15 x 15 mrad^ size 
segmented in three longitudinal sections. The iron return yoke is instrumented with 
streamer tubes read out in projective towers and this provides the hadron calorimetry. 

The tracking and calorimetry information is combined in an energy flow algorithm 
providing a list of objects classified as tracks, photons and neutral hadrons. Particle 
tracks which originate from within a cylinder of 2 cm radius and 20 cm length centred 
on the nominal interaction point, and which have more than 4 TPC hits, are defined 
as good tracks. Hadronic Z decays are selected by requiring the presence of at least 
5 such good tracks and a total visible energy greater than 10% of the centre-of-mass 
energy [|]. 

3 Reconstruction of Charmed Hadron Decays 

Each charm hadron species is reconstructed in the decay mode best suited to a rate 
measurement, i.e. with a large enough and accurately known branching ratio, and a 
sufficiently low background. 




Figure 1: K vr^ invariant mass distribution. The fitted curve is described in section 7. 
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3.1 Meson 



D° mesons are selected through their decay mode D*^ K^tt^. A track with a 
momentum greater than 2.5 GeV/c, to which the kaon mass is assigned, is combined 
with a track of opposite charge with a momentum greater than 1.5 GeV/c, which is 
assumed to be a pion. The contribution due to wrong mass assignment of the two tracks 
from the D*^ is reduced by means of the dE/dx measurement of the kaon candidate 
when available. The dE/dx of the kaon track must be closer to the expectation for a 
kaon than that for a pion: |xk| < IXttI- The pair formed by the two tracks is retained 
if its fractional energy Xe = -E/-Ebeam is greater than 0.5 and if a common vertex with 
a probability greater than 1% is found. The significance of the D° decay length, 
projected on the D° momentum, is required to be greater than 1. The invariant mass 
distribution is shown in Fig. |l]. 
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Figure 2: Distribution of the Kvrvr - Ktt mass difference in four different ranges of the Kyrvr 
energy, after requiring that the Kvr mass is consistent with the mass (histogram with error 
bars). The background fit is also shown. 
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3.2 D*+ Meson 



The D*^-^D°7r+, D°^K^7r+ decay chain is reconstructed by considering triplets of 
charged tracks with total charge ±1. The pion mass is assigned to the two particles 
with the same charge, and the kaon mass to the remaining one. The total energy of 
the triplet is required to be greater than 0.1 x -Ebeam- Any K~7r+ combination whose 
mass IS within 30 MeV/c^ of the D° mass is kept. For such a combination, the mass 
difference M(K7r7r) — M(K7r) is calculated and required to be within 1.6 MeV/c^ of the 
nominal D*+ — D° mass difference. The combinatorial background level is a steeply 
falling function of Xe and background rejection must therefore be enhanced as one 
goes to lower Xe- For this, the scalar nature of the D° and its significant lifetime are 
used as follows. The distribution of the cosine of the D*^ decay angle (defined in the Kvr 
rest frame as the angle 6* of the kaon with the Kvr line of flight in the laboratory) is 
isotropic, whereas the background peaks forward and backward, and so, for Xe < 0.45 
it is required that —0.8 < cos^^* < 0.9. The decay length significance is required 
to be positive for 0.3 < Xe < 0.4 and greater than 1 for Xe less than 0.3. The 
probability for the secondary vertex fit is required to be greater than 1% for Xe < 0.15. 
The mass-difference distributions in four Xe ranges are shown in Fig. 0. 

3.3 D+ Meson 

The meson is reconstructed in the decay mode D+ ^K^7r+7r+. Combinations of 
three tracks are formed with the same mass assignment as for the D*+ reconstruction. 
To fight the high background of this decay, the kaon track candidate is required to have 
a momentum greater than 2.5 GeV/c and an associated dE/dx satisfying Xk + Xtt < 1, 
one of the two pion tracks must have a momentum greater than 1.5 GeV/c, while the 
other must have more than 0.75 GeV/c. Both of the possible Kvr combinations are 
required to have a KTrvr-Kyr mass difference larger than 0.15 GeV/c^, in order to reject 
candidates consistent with the decay D*"*" — > tt'^D'^ tt+K^tt+X. Triplets of tracks 
with a total energy greater than 0.5 x -Ebeam and forming a common vertex with a 
probability greater than 1% are retained. A projected decay length significance greater 
than 1.5 is required for the reconstructed D"*" vertices. The invariant mass distribution 
is shown in Fig. |^. 

3.4 D+ Meson 

The D+ is reconstructed in the decay chain D+ — ^ ^tt"*" with a subsequent decay 
— > K"'"K~. The kaon momenta must exceed 1.5 GeV/c and the pion momentum 
must be greater than 2.5 GeV/c. The kaon candidates are required to have dE/dx 
information available and to fulfill xk + Xtt < 1- In addition all tracks for which 
dE/dx is available must satisfy \xx\ < 2.5 for the given particle hypothesis. If 
|Mkk — M(f)\ < 5 GeV/c^ and > 4 GeV/c the candidate is kept and the three 
particles fitted to a common vertex. Candidates with a greater than 15 (there are 
3 degrees of freedom) are rejected. A cut on |cosA*| > 0.4, where A* is the angle 
between the kaon and the pion in the rest frame, motivated by the P-wave nature 
of the decay, strongly suppresses the background. The mass distribution of the D+ 
candidates is shown in Fig. ^. 
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Figure 4: K+K 7r+ invariant mass distribution, with the fit result. The peak at low mass 
comes from the D"^ 
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3.5 D*+ Meson 



The D*+ meson is reconstructed in the decay channel D*"*" — > D+7, followed by 
D+ K+K~7r+ through either the 07r+ or the K*°K+ decay mode. 

In the D+ reconstruction, tracks are taken as kaon candidates when the dE/dx 
measurement satisfies |xk| < IXttI, and as pion candidates when |xk| > IXttI- To select 
the mesons decaying into 07r+, pairs of opposite-charge kaon candidates satisfying 
|Mkk — Mfp\ < 8 MeV/c^ are required to form a vertex with a pion candidate with 
a probability greater than 1% and | cosA*| greater than 0.4. candidates decaying 
into K*K are selected when kaon and pion candidates with opposite charges fulfill 
\MKn - Mk* \ < 25 MeV/c^ and | cos A*| > 0.6. 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

M(D j)-M(D^) (GeV/c^) M(D j)-M(D^) (GoV/c) 

Figure 5: M(D*~^) — M(D^) distributions with conversions in the bb-enriched sample (a) 
and in the cc-enriched sample (b), and with calorimetric photons (c), (d). The binnings have 
been chosen to match the different resolutions in the two photon detection methods. The 
dashed hne is the background fit described in Section 6 

Photon candidates for the D*+ are selected when they have an energy greater than 
0.6 GeV in the electromagnetic calorimeter and do not form a 7r° candidate. Photon 
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conversions in the detector are identified by searching for electron tracks that are 
geometrically consistent with the conversion hypothesis. To increase the efficiency 
for converted photons, electron-positron pairs where one track is not detected are also 
identified, as described in 0. Photon conversion candidates must have an energy larger 
than 0.3 GeV to be taken into account. 

To be able to measure separately the production rate of D*"*" mesons in cc and 
bb events, two different subsamples of D*+ candidates are selected. One subsample 
is enriched in bb events and consists of D*"*" candidates with a momentum 
8 GeV/c < < 20 GeV/c, a kaon momentum px > 1-5 GeV/c, a pion momentum 
Pt, > 1 GeV/c, a fitted decay length 6i > 500 /xm, and a photon energy < 3 GeV. 
The bb purity of this D*"*" signal is 85%. In the other subsample, enriched in cc events, 
the requirements are p^^ > 20 GeV/c, pk > 2.5 GeV/c, p^, > 1.5 GeV/c, 6i > 0, and 

< 5 GeV, corresponding to a purity of the D*+ signal of about 60% in cc events. 
In Fig. ^the M(D*) — M(Ds) distributions with calorimetric photons and conversions 
are plotted. 

3.6 A+ Baryon 

The A+ is observed via the decay mode A+ — >-pK~7r+. To form a A+ candidate three 
tracks are required, each having at least one vertex detector hit. The momenta must 
be greater than 5.0 GeV/c, 3.0 GeV/c and 1.0 GeV/c for the proton, kaon and pion 
candidates respectively and the total energy must exceed 0.2 x i^boam- For particle 
identification, \xh\ < 2.5 is required for the given particle hypothesis and two pion veto 
cuts are applied: Xtt < —2.0 for the proton candidate and Xtt < —1-0 for the kaon 
candidate. 

After a fit of the tracks to a common vertex the x^ probability has to be greater 
than 1%. The mass distribution of the A+ candidates is shown in Fig. ^ 




Figure 6: pK tt~^ invariant mass distribution with the fit result. 
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4 D*+ Spectrum and Charm Fragmentation 

Because of its large signal/background ratio, the D*+ sample, selected as described in 
section 3.2, is used to measure the fractional energy spectrum. The Xe range between 
0.1 and 1.0 is divided into 18 equal bins. In each bin, an analytical form is fitted 
outside the D*+ peak to the mass- difference distribution, so that the background can 
be subtracted bin by bin. The efficiency-corrected Xe distribution is shown in Fig. |^ 
and the numerical values of the production rates are given in Table |l|. 



Table 1: Differential D*^ production rate as a function of the fractional energy Xe- The 
errors are statistical only. 



Xe range 


1 /N(Z-^hadrons) x (dN(D*±) /AXe) 
(10-3) 




0.10-0.15 


7.47 ±0.63 


0.15-0.20 


9.03 ± 0.49 


0.20-0.25 


10.42 ± 0.44 


0.25-0.30 


10.76 ± 0.43 


0.30-0.35 


9.89 ±0.38 


0.35-0.40 


8.97 ±0.35 


0.40-0.45 


8.17 ±0.32 


0.45-0.50 


6.94 ±0.28 


0.50-0.55 


6.73 ±0.27 


0.55-0.60 


5.56 ±0.24 


0.60-0.65 


4.94 ±0.22 


0.65-0.70 


3.49 ±0.18 


0.70-0.75 


3.13 ±0.17 


0.75-0.80 


2.00 ±0.14 


0.80-0.85 


1.27 ±0.11 


0.85-0.90 


0.50 ±0.07 


0.90-0.95 


0.27 ±0.05 


0.95-1.00 


0.06 ±0.03 



The D*"*" spectrum is interpreted as the sum of three contributions: hadronisation 
of a charm quark produced in a Z — >■ cc decay, decay of a bottom hadron from a Z — bb 
decay, and gluon splitting into a pair of heavy quarks which hadronise or decay into a 
D*"*". The Xe distribution of the three contributions are taken from the Monte Carlo 
simulation which uses the Peterson et al. fragmentation scheme |p for heavy quarks. 
The bb contribution is compared with data using a high purity b-tag in the opposite 
hemisphere (Fig. ^). A possible contamination from two-photon production has been 
investigated by a Monte Carlo study normalised to data [|lOl and found to be negligible 
(< 0.25%). 

To describe the production mechanism, three parameters are fitted to the data : 
the fraction fg of D*"*" from gluon splitting, the fraction of bb events in the remaining 
(no gluon splitting) sample, and the Peterson fragmentation parameter e^- The fit is 
carried out in two steps iteratively. Firstly a two-dimensional fit to the distribution 
{XE.Btag), where ;Btag is the tagging variable described in Section 2 is performed in 
order to constrain fb- Xe is required to be larger than 0.25 to reduce to a negligible 
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level the gluon splitting contribution. The shape of the c and b contributions are taken 
from a full Monte Carlo simulation. The small gluon splitting contribution is taken 
in size and shape from the Monte Carlo simulation. The combinatorial background 
distribution is taken from the side band of the Knir — Kvr mass difference in the data 
sample. The {XE,Btag) plane is segmented into 62 bins of variable size in such a way 
that they all contain enough statistics for a fit of /b to be performed. 
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Figure 7: Distribution of Xe for the D*^'s. The data (points with error bars) are compared 
to the fit described in the text. The three contributions: bb (dotted line), cc (dashed line) 
and gluon spUtting to heavy quarks (dashed-dotted) are shown. 



In a second step, fixing /b to this value, /g and Ec are extracted from a fit to the 
Xe distribution using the function 



+ (1 - Jh) 



A^bb dXE 



Nc-c dXE 



1 dN^^'"'^ 
Ng^cc dXE 



(1) 



The two steps are iterated and once convergence is reached, a minimal of the 2D 
fit of 71.3 for 61 degrees of freedom is obtained for ft, = 0.534 ±0.016, corresponding to 
Rb f{b —>■ D*^)/Rc f{c — > D*+) = 1.15 ±0.06. The uncertainty of this measurement is 
statistically dominated. The fitted distribution of the Btag variable is shown in Fig. ^, 
compared with data. 

The result of the ID-fit is shown in Fig. ^ it has a of 11.9 for 14 degrees of 
freedom (the last three bins are grouped together in the fit). The fitted value of /g is 
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related to the average number of gluon splitting events to charm per Z hadronic decay 
and gives rig^cc = (4.7 ± 1.0(stat.))%. The systematic error on this result however is 
large, due to the efficiency uncertainties at low Xe- In the fit, the JETSET [g| Monte 
Carlo shape for the gluon splitting and the Peterson parameterisation for the heavy 
quark fragmentation functions |]^ are assumed. In the next section, a more sensitive 
analysis dedicated to this measurement is described. 
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Figure 8: Distribution of Xe for the D*"'"'s from the bb source compared with the Monte 
Carlo prediction. The data (points with error bars) are obtained by applying the Ufetime- 
mass btag on the hemisphere opposite to the D*^, and corrected for efficiency. The remaining 
charm contribution is subtracted. 

The fitted value Sc = (33.9 ± 3.7) x 10~^ can be converted into a measurement of 
the average fractional energy of the D*+ in Z — cc events of 

{XeCD*))^- = 0.4878 ± 0.0046(stat.) ± 0.0061(syst.). 

The quoted systematic uncertainty is dominated by the choice of parameterisation 
of the fragmentation function, as in Ref. 

Using the Monte Carlo shape for the very small extrapolation to the entire Xe 
range (from 0.1 to the threshold at 0.044), the number of D*^ per hadronic Z decay in 
the mode considered is found to be 



nD.±x5(D*+ D°7r+)x5(D° K"7r+) = (5.114±0.067(stat.)±0.072(syst.)) x 10 



-3 



10 



A 0.5% uncertainty from the finite Monte Carlo statistics is included in the statistical 
uncertainty. This result includes the gluon splitting contribution. Subtracting this 
contribution and using the ratio of the bb to the cc contribution quoted above, the 
probability for a c-quark to hadronise into a D*+ is measured to be 

/(c ^ D*+) = 0.2333 ± 0.0102(stat.) ± 0.0084(syst.) ± 0.0074(B.R.). 

The normalisation of this result relies on the branching fraction given in ref. 0]. 




Figure 9: Distribution of the b-tagging variable for the D*='='s with Xe > 0.25 The data 
(points with error bars) are compared with the projection onto the Btag axis of the fit 
described in the text. The two contributions: bb (dotted line) and cc (dashed line) are 
shown. 

Systematic errors on these quantities arise from various sources. The uncertainty 
from /b has been evaluated by varying it by ±1 standard deviation. The uncertainty 
on the efficiency has been estimated on a cut by cut basis. The uncertainty from the 
mass and mass-difference cuts have been assessed by varying the fitted mass and mass 
difference within errors yielding a 0.5% and 0.3% contribution to the relative error on 
the efficiency. The Prob(x^) cut efficiency has been compared between data at high 
Xe and Monte Carlo, and found to differ by 8.3%. This is taken as the uncertainty on 
the efficiency of this cut and results in a 0.6% uncertainty on the overall normalisation. 
To assess the uncertainty from the cut on decay length significance, the resolution on 
this quantity has been varied by 10%. The resulting relative change in efficiency is 
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found to be (1.6 ± 0.6)% for the 51 /a > cut and (0.06 ± 0.90)% for the 6i/a > 1 cut. 
Taking a 1.6% uncertainty in the range 0.3 < Xe < 0.4 and 0.9% for Xe below 0.3, the 
resulting uncertainty on the overall normalisation is 0.4%. The probability of nuclear 
interaction in the detector is simulated at the 10% level, determined from a study of 
tau pairs where the pion from one of the taus undergoes such an interaction in the 
detector. This yields a 0.7% relative error on the efficiency. For the systematics from 
background subtraction, the analytical shape of the background has been modified and 
the fit range for the mass-difference fit has been varied. This results in a 0.2% overall 
error on the absolute normalisation, localised in the first three Xe bins where the 
background is large (the uncertainty reaches 4.5% in the first bin). Summing linearly 
the errors from the P(x^) cut and the S£/a cut, as they both pertain to vertexing 
and are strongly correlated, and then summing in quadrature all the contributions, the 
overall normalisation uncertainty is 1.4%. Details are given in Table |^. As a check, 
the whole analysis has been repeated with the decay length significance cut relaxed by 
1 unit and the cut in probability for the D° vertex removed. The final number of 
events changes by 1%. 



Figure 10: Distribution of Xe for the D* s from the cc source. The data (points with error 
bars) are compared to the Monte Carlo simulation. 

Some of these systematic effects are Xe dependent. They have been propagated to 
the gluon splitting measurement and the average Xe measurement. 

As a check on the stability of the result, the fit has been repeated with /b left free. 
The fitted parameters are perfectly consistent with the result of the whole procedure. 
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but the uncertainty on is increased by a factor of 1.7. 

In Fig. |l^, the spectrum of the D*s from charm is shown. This spectrum (points 
with error bars) is obtained from the data points of Fig. |^ or Table |I], from which the b 
contribution and the gluon sphtting contribution are subtracted. The b contribution is 
extracted from the b-tagged data of Fig. |, normahsed using the fitted /b. The Monte 
Carlo shape, normalised to the measurement presented in the next section, is used for 
the gluon-splitting contribution. 

Table 2: Relative systematic errors in percent on the D*^ production rate measurement. 
Some of them are strongly Xe dependent (see text). The two vertexing uncertainties are 
considered as fully correlated and added linearly. 



Source 


5e/e 


AM cut 


0.3% 


M(K7r) cut 
P(x') cut 


0.5% 


0.6% 


6i/cr cut 


0.4% 


nuclear interactions 


0.7% 


Backg. subtraction 


SN/N = 0.2%. 


Total systematics 


1.4% 



5 Measurement of the Gluon Splitting Rate to cc 

In order to determine more precisely the rate Ug^cc one more variable, related to 
the event topology, has been used. In a g — ^ cc event the two quarks from the 
gluon splitting, together with one of the primary quarks, tend to end up in the same 
hemisphere, resulting in a large invariant mass for this hemisphere and a relatively light 
opposite hemisphere. Denoting by Mheavy (^light) the invariant mass of the heavier 
(lighter) hemisphere in an event, the difference AM^ = Mheavy — ^light has been used 
to discriminate between D*'s produced from primary quarks and from gluons. The 
distributions of AM^, as given by the JETSET event generator are shown in 
Fig. 11. In JETSET D*'s from primary quarks are nearly evenly shared among the 
hemispheres while ~ 75% of the D*'s from gluons are found in the heavy hemisphere. 

The sample of reconstructed D*"*" mesons is divided into two samples depending 
upon whether the D*"*" is in the light or heavy hemisphere. A two-dimensional binned 
likelihood fit is then performed to the distributions in the Xe/ AM^ plane of these two 
samples using the function : 



(1 - h) 



/b 



bb 



1 c/2jYg^cc 



1 



N^-c dXEdAMu 



iVg^ec dXEdAMu 



(2) 



For this analysis the event selection has been tightened in order to ensure that the 
event is fully contained in the detector: The total energy Ech of all charged tracks 
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is required to exceed 15 GeV and events strongly unbalanced in momentum projected 
along the beam axis are rejected by requiring | X^PzlZ-^ch < 40%. This selection has 
an efficiency of 89.3% for hadronic Z decays without the process g — > QQ and 91.9% 
for g — i> QQ events (Q stands for b or c), as estimated from Monte Carlo. 
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Figure 11: Difference of heavy and light hemisphere masses in Z decays with (dashed line) 
and without (full line) gluon splitting to heavy quarks (from JETSET). The events entering 
these distributions are required to contain at least one D*^. Q stands for b and c, and q for 
all flavours. 

The D*"^ mesons are reconstructed with cuts similar to the ones used in section 3. 
The hemispheres are defined by a plane perpendicular to the thrust axis and their 
invariant masses are calculated using energy flow objects with an energy above 1 GeV. 
Again, the analysis proceeds in two steps. A fit to the one-dimensional Stag distribution 
yields the fraction fstag = 0.524 ± 0.015 of D* from b hadron decays in the selected 
sample. The error includes the systematic uncertainties but is dominated by the 
statistical error. Then the two two-dimensional distributions Xe/ ^M^i depending 
upon whether the D* is in the heavy or light hemisphere, are fitted simultaneously 
using a binned likelihood fit (18 x 18 bins). The b-contribution excluding gluon splitting 
/b(l — /g) is constrained to fstag by a term. The result of the fit is 

/g = 0.0522 ± 0.0077 , /b = 0.5308 ± 0.0098 



from which n„^cc is computed as 



n. 



g^cc 



Rr 



^g-^cc f 



/b 



1 e„ 



■>bb 



1 ~ /g Rh Cb 



+bb 



(3) 
(4) 



= (3.23±0.48(stat.) 
Here, is the efficiency to reconstruct a D*+ from a given source i. The term in brackets 



provides a ~ 10% correction for the g ^ bb contribution which is taken from ref. |11]: 
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Figure 12: Hemisphere mass difference distributions. The data (points) axe subdivided into 
four bins with D* in the heavy/hght hemisphere and for low/high Xe- The fitted contribution 
from primary quarks is shown as the soUd histogram, the g QQ contribution is added as 
cross-hatched area. 
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?T-g^bb = (2.77 ± 0.42 ± 0.57) x 10^'^. In Fig. |T2| projections of the two-dimensional 
distributions on AM^ are shown for the two hemispheres for two ranges in Xe- The 
agreement between data and Monte Carlo is good, in particular it can be seen that the 
high-X^; data are well described without the gluon splitting contribution but that its 
inclusion is necessary for Xe < 0.35 in the heavy hemisphere. The probability for 
this distribution degrades from 0.40 to 0.001 for the four bins when the gluon splitting 
contribution is neglected. In order to show more clearly the shape of the gluon splitting 
signal in this region, Fig .|13] shows the data with the fitted contribution from primary 
quarks subtracted. Good agreement with the Monte Carlo prediction is observed. 

Systematic errors have been studied for the efficiency of D* reconstruction, the 
fitting procedure, the modelling of the Monte Carlo distributions, and the parameters 
entering Eq. (|^) . They are summarised in Table |^. 
The errors on the number of reconstructed D*"'"'s have been estimated as in the 



Table 3: List of absolute systematic errors on n. 



Source 


Ang^ec (10-2) 


D*^ reconstruction 


0.25 


Monte Carlo statistics 


0.19 


b-contribution 


0.09 


Hemisphere selection 


0.15 


£c 


0.15 


£b 


0.19 


XE{h D*) from data 


0.15 


D** ^D*X 


0.04 




0.15 


Hard gluon 


0.08 


Herwig 


0.11 


Ariadne 


0.21 


Total 


0.53 



previous section by comparing the efficiencies of the cuts between Monte Carlo and 
data and varying them within their uncertainties. The normalization of the background 
has also been changed within the errors. 

The finite Monte Carlo statistics has been accounted for by repeating the fit with 
the distributions randomly smeared. The fit has also been redone not applying the 
constraint on the b contribution, and leaving free the relative fractions of D*'s in the 
two hemispheres for the c and b contributions separately. 

Several parameters affecting the shape of the Monte Carlo distributions have been 
varied: has been changed within the errors derived from the above measurement 
and within eb = 0.0045 ± 0.0009. The Xe shape for b — D* as derived from the 
b-tagged sample has been used for the fit. The dependence of the fraction of D* mesons 
originating from a D** has been estimated by changing this rate by ±30% for b and c 
separately. 

The effects of gluon emission in the parton shower have been assessed by changing 
the JETSET parameter Aps such that it corresponds to a ±4% variation in ^^(Mz). 
The rate of events in which both primary quarks recoil against a hard gluon and 
populate the same hemisphere has been varied by 30%. Finally, the AMh shapes as 
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given by HERWIG 5.8 and ARIADNE 4.08 g have been used. For ARIADNE, 
the transverse momentum of a splitting gluon was required to be larger than its 
virtuahty, as suggested in Ref. []TB[. ARIADNE yields the largest discrepancy and 
this has been taken as the systematic error on the Monte Carlo model. Adding all 
systematic errors quadratically, the average fraction of hadronic Z decays where a 
gluon splits in a cc pair is found to be 



n, 



g^cc 



(3.23 ± 0.48(stat.) ± 0.53(syst. 
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Figure 13: Comparison of the data with the fitted contribution from primary quarks 
subtracted (points) and the predicted gluon splitting distribution (histogram). 



6 Measurement of the T)*^ production rate 



The selection efficiencies for finding D*^ mesons are evaluated from the Monte Carlo 
using the mix of D^ — > 07r^ and D*+ — > K*^K~^ decays predicted by the branching 
ratios in [|14|. They are given in Table |^ for the selection described in Section 3.5. 



Table 4: Selection efficiencies for D*^ mesons (the errors are statistical only). 





conversions 


calorimetric photons 




Z ^ cc (%) 


Z ^ bb (%) 


Z ^ cc (%) 


Z ^ bb (%) 


cc-enriched 


0.38 ±0.04 


0.12 ±0.02 


2.72 ±0.12 


0.85 ±0.07 


bb-enriched 


0.10 ±0.02 


0.29 ±0.03 


0.57 ±0.05 


1.60 ±0.09 



The number of events in the data samples was evaluated by fitting the D*"*" spectrum 
with the sum of two Gaussian functions for the signal and a polynomial for the 
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background. In the case of the calorimetric photons, the means and widths of the 
two Gaussians and the ratio of their area are taken from Monte Carlo. For converted 
photons the means of the two Gaussians and the width of the broad one are left free 
in the fit. The background is parameterised to reproduce the Monte Carlo. From 
the fit the number of events in the cc-enriched sample is iVc = 26 ± 7 (139 ± 22) and 
in the 66-enriched is A^^f, = 36 ± 11 (106 ± 37) with converted (calorimetric) photons. 
The production rates /(c Dg^) and f{b Dg^) can be calculated by solving two 
equations of the following form, corresponding to the cc and bb enriched samples: 

iVD*+ = 2iVqq X i?(D:+ ^ D3+7)5(D3+ ^ 07r+)5(0 ^ K+K-) 

X [RJ(c ^ D:+)e, + Rb/(b D:±)eb], (5) 

where A^'qq is the number of hadronic events in the data sample, ec (eb) the selection 
efficiency for D* mesons in cc (bb) events and /(c — > D*) (/(b D*)) the probabilities 
of a c (b) quark giving rise to a D* meson. A similar formula can be written for the 
— > K*''K"'" channel. The results from the whole data sample are: 

/(c ^ D*+) = 6.9 ± 1.8 (stat.) ± 0.7 (syst.) ± 1.8(BR)% 
/(b ^ D*+) = 11.3 ± 3.5 (stat.) ± 1.0 (syst.) ± 2.8(BR)% 

Here the last error is due to the final state branching ratios taken from [|I4|, while the 
systematic error includes the uncertainty in the selection efficiencies and background 
parameterisation in the fit. The latter contribution is evaluated using different 
functional shapes for the background; a further check is made by taking the background 
shape evaluated bin by bin from Monte Carlo and using it in the fit. 

7 Production Rates of the Ground States 
7.1 Production of D° and D+ 

The number of D° and present in the samples are extracted from a fit to the 
invariant mass distributions (Figs |l|, ^), where the signal is parameterised by two 
gaussians with a common mean and the combinatorial background by a polynomial 
function. Resonant background contributions, such as D° — >■ K~K"'", — >■ ^vr"*", 
D+ K*K, and D°-^K-7r+ where the two mass assignments are reversed, are taken 
into account in the fit. The fractions of D° and D"*" originating from the process 
Z ^ cc are measured by applying, on the hemisphere opposite to the reconstructed D 
meson, the lifetime-mass tag selecting b hemispheres with 99% purity. The charm 
fraction in the D meson samples is then given by fc = where VTa-tag is the 

fraction of D mesons surviving the b-tag cut. The b hemisphere b-tag efficiency eeb 
is evaluated on data using a double-tag technique and esc, the c hemisphere b-tag 
efficiency, is obtained from Monte Carlo simulation ||^. Both efficiencies are corrected 
for the selection bias induced by requiring the presence of a high momentum D meson 
in the opposite hemisphere. The fractions of D° and mesons from charm are found 
to be: /c(DO) = 0.787 ± 0.019 ; /c(D+) = 0.797 ± 0.020, where the quoted errors are 
statistical. 

The small fraction of D mesons from gluon splitting has been deduced using the 
measurement from section 5, the ALEPH measurement ?T-g^bb = (2.77 ±0.42 ±0.57) x 
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10 ^ pJTl, and the efficiency ratio for the selection of a D meson obtained from 
simulation. 

The number of fitted D mesons, together with the charm fractions and the gluon 
sphtting contamination, are reported in Table |^. 



Table 5: D'^ and selection: number of candidates; fraction of c — > D mesons in the 
sample; efficiency of c ^ D meson reconstruction; percentage of D mesons from gluon 
splitting. All quoted errors are statistical. 





N. candidates 


Charm fraction 


Efficiency 


QQ(%) 




7871 ± 129 


0.787 ±0.019 


0.200 ±0.004 


0.8 


D+ 


7409 ± 153 


0.797 ±0.020 


0.191 ±0.002 


0.8 



From the numbers in Table ^, the following branching fractions are obtained: 
Rc X /(c ^ D°) X fi(D° ^ K"7r+) = (0.370 ± 0.011) x 10"^ 

Rc X /(c ^ D+) X 5(D+ K-7r+7r+) = (0.368 ± 0.012) x 10"^ 

Different sources of systematic uncertainty have been investigated and are given in 
Table M. The ffist three errors are related to the detector performance. The fraction of 



Table 6: Relative systematic errors in percent on the D*^, D+, Dg^ and production rate 
measurements . 



Source 


DO (%) 


D+ (%) 


Ds+ (%) 


A+ (%) 


Nuclear interact. 


0.5 


0.7 


0.3 


0.3 


dE/dx 


0.9 


0.4 


1.2 


2.2 


vertexing 


3.8 


2.3 


1.7 


1.7 


charm fraction 


0.7 


0.7 






gluon splitting 


0.3 


0.3 


0.3 


0.3 


fragmentation 


3.7 


3.7 


1.9 


4.0 


charm physics 


1.6 


2.1 






ffiting funct. 


1.5 


1.5 


4.0 


1.6 


MC Stat. 


2.0 


1.0 


3.0 


2.0 


Total systematics 


6.2 


5.3 


5.8 


5.5 



tracks with dE/dx information has been compared between data and Monte Carlo with 
a sample of candidates selected in the side-band regions of the D° and D+. The Monte 
Carlo has been corrected for differences and the uncertainty on the ratio data/MC 
has been taken as a systematic error. The dE/dx calibration and resolution have been 
studied on data for a pure kaon sample from the decay mode D*"*" — D^vr"*", D'' — K~7r+. 
Data/MC differences have been taken into account to obtain the dE/dx cut efficiency 
and the associated systematic error. The fraction of D mesons forming a common 
vertex with a x^-probability greater than 1% is obtained from data, from the fraction 
of selected candidates under the D meson mass peak, corrected by the fraction of 
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background candidates estimated on sideband regions. The resolution on the decay- 
length measurement has been checked on D° and D+ sideband samples, where most 
of the candidates carry no lifetime information. The difference found between data 
and Monte Carlo has been taken as a systematic error for the decay length significance 
cut together with the error coming from the D'' and lifetime uncertainties. The 
fraction of tracks not reconstructed due to nuclear interaction in the D° — > K~7r"*" and 
D"*" — > K~7r+7r+ decays is estimated from Monte Carlo to be 5 % and 7 % respectively, 
with a systematic error of 10 %. 

The main systematic error on the determination of the charm fraction comes from 
the correction factor applied to the b-tag efficiencies. This arises from the selection of 
an energetic D meson in one hemisphere which favours events with no gluon radiation, 
hence increasing the efficiency for the b-tag in the opposite hemisphere. A comparison 
between data and Monte Carlo of the spectrum of the most energetic jet in the 
hemisphere opposite to the D meson has been performed to check the simulation of 
gluon emission. The quoted error on the correction factors arises from the limited 
Monte Carlo statistics and the difference between data and Monte Carlo for gluon 
emission; this gives a relative error on the charm fraction of 0.7%. The relative 
systematic error due to the uncertainty on the b-tag efficiencies is very small, of the 
order of 0.2%. The systematics related to the charm fragmentation and the simulation 
of the decay chain leading to a D° and D+ are evaluated by varying the fragmentation 
function and the fraction of charmed vector, pseudoscalar and excited meson (10 to 
30%) states produced in the charm hadronisation. The Peterson fragmentation function 
was varied to reproduce the measured D*"*" Xe spectrum for different fractions of 
produced D meson states. Different kinds of fit have been performed to assess the 
systematic error coming from the fitting. For the signal function, the mean and width 
of the narrow gaussian distribution is left free, while the relative normalisation and 
width of the second gaussian is fixed to the Monte Carlo expectation. The fixed 
parameters have been varied within 5% and the fit repeated with other background 
paramctcrisations (second-order polynomial, exponential). Fitting the background 
outside the mass peak and counting the number of candidates above the mass peak 
gives a consistent result. 

The final result on the D° and D+ production rate measurements in Z — > cc is 

Rc X /(c D°) X S(D° ^ K-7r+) = (0.370 ± 0.011 ± 0.023) x 10"^ 

Rc X /(c D+) X fi(D+ K-7r+7r+) = (0.368 ± 0.012 ± 0.020) x 10"^ 

Dividing these numbers by the Standard Model value of Rc and the branching 
fraction for the reconstruction mode, the following fractions are obtained: 

/(c ^ D°) = 0.559 ± 0.017 ± 0.0035 ± 0.013(BR) 

and 

/(c ^ D+) = 0.2379 ± 0.0077 ± 0.0129 ± 0.0190(BR). 

The consistency between the rate measurements is investigated. Under the 
assumption that the observed difference between /(c ^ D°) and /(c ^ D+) is only 
due to D* production, i.e. if primary D's and D*'s are produced evenly in the two 
charge states, the following ratio must equal unity: 

/(c^D°)-/(c^D+) 
2 /(c ^ D*+) x B{D*+ 0%+) 
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This ratio is measured to be 1.02 ± 0.12. 



7.2 Measurement of the Py Ratio 

The ratio Py = where V is the fraction of D* (vector) produced and P the 

fraction of D (pseudoscalar) , after all decays of heavier excited states, can be derived 
either from the ratio /(c D'^)//(c D^), or from the D*^ and rates. The latter 
leads to 

Py = 0.595 ± 0.045. 

This value is more consistent with the predictions of the thermodynamical 
approach [|I^ and the string fragmentation approach of ref. |T^, which both predict 



0.66 for Py, than from the naive expectation of 0.75 from counting the spin states. 
From the analysis of section 6, the Py ratio is found to be 

0.60 ±0.19. 



for and D*+ mesons. 



7.3 Study of A+ and D+ production 

The c fractions in the and A'^ samples are extracted by the same method as for 
the D° and D+. The background fit includes the Monte-Carlo shape and the proper 
normalisation of contributions from D° and D"^ to the invariant mass distributions. The 
efficiency correction is similar, except that the data are counted in 8 Xe bins between 
0.2 and 1, allowing the fit of the Xe distributions to a Peterson-based shape for the 
purpose of extrapolation down to the threshold. The breakdown of the systematic 
errors is very similar to the one of the previous section and is given in the last two 
columns of table As the relative bb to cc contributions are fitted to the data, the 
uncertainty on the charm fraction is included in the statistical error. 

The following and product branching ratios in Z — > cc are obtained: 

Rc X /(c ^ D3^) X 5(D+ ^ (/)7r+) x 5(0 ^ K+K") = (0.352 ± 0.057 ± 0.021) x 10"^ 

Rc X /(c ^ A+) X 5(A+ ^ pKvr) = (0.673 ± 0.070 ± 0.037) x 10^1 

Dividing these numbers by the Standard Model value of Rc and the branching 
fraction for the reconstruction mode, the following fractions are obtained: 

/(c ^ D+) = 0.116 ± 0.019 ± 0.007 ± 0.030(BR) 

and 

/(c A+) = 0.079 ± 0.008 ± 0.004 ± 0.020(BR). 

Alternatively, one can use the sum of Rc x /(c — > Xc) over all the fundamental 
states of the various hadron species to determine Rc as in section 8. 
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8 Measurement of Rc from Charm Counting 



By dividing the values of the product branching fractions found in previous sections by 
the corresponding decay branching fraction, taken from the following individual 
charm meson rates are obtained. 

Rc X /(c ^ D°) = 0.0961 ± 0.0029 ± 0.0060 ± 0.0023(5i?) 

Rc X /(c ^ D+) = 0.0409 ± 0.0013 ± 0.0022 ± 0.0033(5/2) 

Rc X /(c ^ D+) = 0.0199 ± 0.0032 ± 0.0012 ± 0.0050(5/?) 

Rc X /(c ^ A+) = 0.0135 ± 0.0014 ± 0.0007 ± 0.0035(5/?) 

The small additional contribution from Sc and Qc is estimated to be 0.0034 with 
a 50% uncertainty using usual stangeness suppression factors. Summing all these 
contributions, a value of 

Rc = 0.1738 ± 0.0047(stat.) ± 0.0088(syst.) ± 0.0075(BR) 

is found where the systematic uncertainty is a linear sum of the common 
contributions (nuclear interactions, dE/dx, vertexing, gluon splitting, charm fraction 
and fragmentation), added quadratically to the other contributions. This is in good 
agreement with the Standard- Model expectation 0.1719 |^ and the accuracy (7%) of 



this measurement is similar to others ^ . The breakdown of the systematic errors 
on Rc is given in Table 0. 



Table 7: Systematic errors on Rc. 



Source 


5Rc 


MC statistics 


0.00207 


Charm physics 


0.00176 


Fragmentation 


0.00599 


Charm fraction 


0.00096 


Gluon splitting 


0.00051 


dE/dx 


0.00156 


Nuclear interactions 


0.00087 


Vertexing 


0.00516 


Mass fits 


0.00177 


5(D0 ^ K-71+) 


0.0023 


5(D+ ^ K-7r+7r+) 


0.0033 


5(D+ ^ 07r+) 


0.0050 


5(A+ pK~7r+) 


0.0035 


baryons not decaying 




to A+ 


0.0017 


Total internal 


0.0088 


Total external 


0.0075 
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9 Conclusions 



A significant improvement of tlie knowledge of charm production in Z decays has been 
achieved. This is due, with respect to previous analyses, to a better rejection of the 
combinatorial background and a better b/c separation, both coming from the use of 
the vertex detector. The D*^ spectrum is measured and interpreted as a sum of three 
contributions: hadronisation of a charm quark produced in a Z — cc decay, decay of a 
bottom hadron from a Z — bb decay, and gluon splitting into a pair of heavy quarks 
which hadronise or decay into a D*^. This scheme is found to describe the data very 
well. 

The probability for a c quark to hadronise into a D*+ meson is found to be 

/(c ^ D*+) = 0.233 ± 0.010(stat.) ± 0.011(syst.). 

A fit to the spectrum yields the average fractional energy of the D*^ in Z ^ cc 
events 

{XeCD*))^- = 0.4878 ± 0.0046 ± 0.0061. 

Using the heavy - light hemisphere mass difference distribution the average number 
of gluon splitting events to charm quark pairs per Z hadronic event is found to be 

ng^cc = (3.23 ± 0.48 ± 0.53)%. 

The production rates of the main charmed ground states have been measured. The 
effective ratio of Vector/ (Vector+Pseudoscalar) production rates in charmed mesons 
is found to be 0.595 ± 0.045, confirming the rather low value found earlier 
but consistent with recent semi-phenomenological models. For the charmed-strange 
mesons, a value of Py = 0.60 ± 0.19 is found. 

Summing the contributions of all the fundamental charmed states, and including a 
contribution from baryons not decaying to A+, a measurement of Rc is achieved: 

Rc = 0.1738 ± 0.0047(stat.) ± 0.0088(syst.) ± 0.0075(BR). 



This result, combined with other ALEPH measurements |^2|, leads to 

R, = 0.1698 ±0.0069. 
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